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This  report  addresses  two  separate  aspects  of  the  effects  of  the  lower 
ionosphere  on  ELF/VLF/LF  propagation.  First,  full-wave  calculations  are  per- 
formed to  assess  the  effect  of  a ledge  of  ionization  in  the  55-km  to  65-km 
altitude  range  (e.g. , vC-layer^i  on  long-path  propagation.  Second,  relation- 
ships among  wave  frequency,  incidence  angle,  and  refraction  heights  are  derived. 


1.  I have  reviewed  the  Final  Report  by  Pacific  Sierra  Research 
Corporation,  titled  ANALYSIS  OF  ELECTRON  DENSITY-PROFILES  IN 
THE  LOWER  IONOSPHERE.  The  report  describes  work  performed  under 
contract  to  RADC/EEP,  and  represents  the  product  of  that  work. 

2.  The  report  describes  the  effects  of  a low  altitude  ledge  of 
ionization  ('C-layer")  on  the  propagation  of  ELF./VLF/LF  radio 
waves;  and,  derives  relationships  between  wave  frequency, 
incidence  angle,  and  reflection  heights  which  can  be  used  to 
deduce  grazing  incidence  ionospheric  reflectivity  information 
from  steep-incidence  data.  The  research  described  in  the  report 
contributes  to  the  understanding  of  the  complex  interaction 

of  long  radio  waves  in  an  ionized  medium. 

3.  All  aspects  of  the  work  contracted  for  by  RADC/EEP  have  been 
completed  and  described  in  the  subject  report  by  Pacific  Sierra 
Research  Corporation,  and  meets  the  approval  of  this  office. 
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I.  INTRODUCTION 


This  report  comprises  two  independent,  self-contained  parts:  Section 
II  examines  the  effect  of  an  ionospheric  C-layer  on  long-range  ELF/VLF/LF 
propagation;  Sec.  Ill,  the  relationships  among  Incidence  angle,  frequency, 
and  VI.F/LF  reflection  heights.  Subsections  in  both  Sec.  II  and  Sec.  Ill 
present  background  discussions,  derivations,  and  conclusions.  Both  sections 
can,  therefore,  be  treated  as  essentially  separate  reports. 


1 

1 1 • EFFECT  OF  AN  IONOSPHERIC  C-LAYER 
ON  LONG-PATH  Ki.F/VLF/LF  PROPAGATION 

Roughly  speaking,  the  term  "C-layer"  refers  to  Ionization  that 
occurs  below  the  normal  ionospheric  D-laver;  i.e.,  below  altitudes  of 
around  n*<  km.  Many  authors  have  reported  experimental  evidence  for  the 
existence  of  such  a C-layer.  Most  of  this  evidence  is  based  on  steep- 
incidence  ionospheric-ref  lection  data  measured  at  frequencies  between 
IS  kHz  and  100  kHz.  The  characteristics  of  the  C-layer  are  highly 
variable,  depending  on  time  of  day,  season,  location,  and  phase  of  the 
sunspot  cycle.  There  appear  to  be  many  locations  or  times  at  which  a 
detectable  C-layer  does  not  exist. 

As  indicated  above,  most  of  the  radio  data  pertaining  to  the  C-layer 
were  measured  on  transmission  paths  of  no  longer  than  200  km  to  300  km, 
where  the  incidence  angle  at  the  ionosphere  is  not  greater  than,  say  60 
degrees.  For  such  short  paths,  the  propagation  may  be  adequately  repre- 
sented by  the  superposition  of  a ground  wave  plus  first-  and  second-order 
"skv  waves"  that  are  reflected  once  or  twice,  respectively,  from  the 
ionosphere.  At  much  greater  distances,  this  simple  plane-wave  super- 
position becomes  cumbersome,  because  numerous  sky  waves  must  be  retained, 
and  representation  of  the  signal  in  terms  of  a few  waveguide  modes  be- 
comes preferable.  This  section  calculates  the  effect  of  model  C-lavers 
on  Fl.F/VLF/LF  signals  for  pathlengths  so  large  that  the  waveguide-mode 
representation  should  be  used.  We  assume — without  justification — 
that  the  C-layer  exists  over  the  entire  length  of  the  propagation  path. 


IONOSPHERIC  M0DK1.S 


Our  approach  to  assessing  the  sensitivity  of  long-path  Kl.F/Vl.F/l.F 
propagation  to  the  presence  of  a C-layer  consists  of  two  steps.  First, 
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we  make  field-strength  calculations  based  on  a nominal  model  ionosphere 
that  does  not  exhibit  a C-lnyer.  Second,  we  alter  this  model  bv  super- 
imposing sample  C-layers,  recalculating  field  strengths,  and  comparing 
[.  the  results  with  those  obtained  for  the  nominal  model. 

For  our  nominal  model,  we  use  the  DNA  standard  daytime  model  iono- 
sphere ( Knapp  and  Sahuartz , 1975).  Figure  1 shows  the  corresponding 
electron  and  positive  ion-density  height-profiles,  and  N+,  respec- 
tively. Also,  we  assume  the  ion-collision  frequency  to  be  1/40  of  the 
electron  collision-frequency  at  all  altitudes.  The  masses  of  both 
positive  and  negative  ions  are  assumed  to  be  65  AMU  below  40  km,  40  AMI' 
above  65  km,  and  are  linearly  interpolated  between  40  km  and  65  km.  For 
the  model  ionospheres  used  in  this  report,  the  ions  affect  ELF  propaga- 
tion somewhat,  but  VLF  propagation  is  totally  dominated  by  electrons. 

A large  number  of  C-layer  profiles,  measured  by  numerous  investi- 
gators, are  reviewed  and  presented  by  Alpert  (1975),  Bain  (197-1), 
Krasnushkin  and  Federov  (1966),  and  Risbeth  and  Garriot  (1969),  who  also 
give  extensive  bibliographies.  These  experimentally  determined  layers 
exhibit  widely  varying  characteristics,  depending  on  local  conditions 
at  the  time  that  the  data  were  taken.  Moreover,  although  consistent 
with  low-frequency,  steep- Inc idence , reflection  data,  those  layers  are 
not  necessarily  unique;  i.e.,  some  other  assumed  layer  characteristics 
might  also  be  consistent  with  the  reflectivity  data. 

Of  all  the  electron-density  height-profiles  given  in  the  above 
references,  the  most  pronounced  C-layers  had  a maximum  density  of  about 
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100  el/cm  at  a height  of  about  55  km.  Such  a layer  is  .assumed  here, 
and  is  shown  in  Fig.  1 superimposed  on  the  standard  ionosphere.  Through- 
out this  report,  the  C-layer  shown  in  Fig.  1 will  be  called  the  "Bain" 
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C-layer  because  it  corresponds  closely  to  the  mo«f  pronounced  of  the 
models  considered  by  Bain  (1974),  although  a very  similar  profile  was 
also  used  by  Krasnushktn  and  Federov  (1966) . The  C-layer  shown  in 
Fig.  1 is  based  on  measurements  taken  at  summer  noon  in  Great  Britain. 

To  assess  the  effects  of  C-layer  height  on  long-path  propagation, 

3 

we  assume  a ledge  of  density  100  el/cm  and  thickness  6 km  centered  at 
three  heights:  58  km  (Model  A),  55  km  (Model  B) , and  hi  km  (Model  C) . 

Figure  2 shows  the  superposition  of  Model  A on  the  standard  ionosphere. 

Since  Models  B and  C consist  simply  of  moving  the  ledge  of  Model  A down 
or  up,  respectively,  by  3 km,  they  are  not  shown  in  the  figures. 

Although  the  upper  and  lower  boundaries  of  the  C-layer  ledge  are  shown 
as  being  totally  abrupt  in  Fig.  2,  they  are  "smeared"  over  about  1 km 
in  performing  the  numerical  calculations  because  of  the  1-km  height- 
resolution  used  for  the  input  data. 

EQUATIONS  AND  MODE  PARAMETERS 

The  detailed  equations  describing  ELF/VLF/LF  propagation  were  given 
by  Field,  et  al.  (1976)  and  will  not  be  repeated  here.  These  equations 
are  solved  numerically,  accounting  for  the  vertical  inhomogeneity  of  the 
Ionosphere  and  the  curvature  of  the  earth,  thus  obtaining  the  results  given 
below.  The  geomagnetic  field  of  the  earth  is  omitted  from  the  calcula- 
tion, which  is  permissible  for  the  daytime  propagation  conditions  consi- 
dered. To  define  the  notation  and  illustrate  the  key  dependences,  two 
equations  for  the  spatial  dependence  of  the  fields  are  recapitulated 
below . 

We  consider  vertically  oriented  VLF  transmitters  for  which  the  fields 
are  composed  of  a superposition  of  so-called  TM  waveguide  modes.  We  express 
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the  results  In  terms  of  the  magnetic  intensity,  H,  which  is  horizontally* 
oriented  and  is  given  by  Field,  et  al . (7 976): 
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s ind/a 
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*5  -nd/8 . 7x10  -1  n 
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In  Sq . (1),  If  is  the  effective  electric-dipole  moment  of  the  trans- 
mitting antenna;  A,  the  free-space  wavelength;  d,  the  distance  from  the 
transmitter;  a,  the  earth's  radius;  and  c,  the  vacuum  speed  of  l iglit . 

The  quantity  S ^ is  essentially  the  eigenvalue  of  the  a^-waveguide 

mode  and  must  be  computed  numerically.  At  VLF.  however.  S turns  out 

« 

u 

to  have  a magnitude  close  to  unity,  and  the  term  S in  F.q.  (1)  does  not 

u 

appreciably  Influence  the  field. 

Thus,  the  magnitude  of  the  vertical  electric  field  depends  on  the 
state  of  the  ionosphere  through  two  parameters:  A , the  excitation  factor 
for  a vertical  dipole;  n,  the  attenuation  rate  in  decibels  per  megameters 
of  propagation  (dB/Mm) . The  phase  of  the  signal  is  governed  by  the 
relative  phase  velocity,  v/c,  which  must  also  be  determined  numerically. 

Contrary  to  the  situation  at  VLF,  an  ELF  transmitting  antenna  would 
be  horizontally  oriented  and  operated  in  the  cndfire  mode  to  excite  the 
quasi-TM  mode,  which  is  the  only  one  that  propagates  well.  The  appro- 
priate expression  for  the  magnetic  intensity  is 

lltf  JSFiT  A -2n  led/  Av  -qd/ (8 . 7 x 1 (/’) 

1 C A/m  , (2) 

sin  d/a 


More  precisely,  we  consider  the  azimuthal  component,  H where  if 
is  the  usual  azimuthal  coordinate  in  polar  spherical  coordinates. 


(1) 


8 


whore  f is  the  wave  frequency,  0 Is  the  effective  ground  conductivity 

in  the  vicinity  of  the  transmitter.  A,  is  the  excitation  factor  for  a 

h 

horizontal  dipole  in  the  endfire  direction.  As  for  VLF,  the  mode  para- 
meters to  be  calculated  are  the  excitation  factor,  A.  , the  attenuation 

n 

rate,  n.  and  the  relative  phase  velocity,  v/c. 

As  defined  in  Eqs.  (1)  and  (2),  A and  A,  have  units  of  inverse 
distance.  In  the  idealized  limit  of  a perfectly  reflecting,  sharply 
hounded  ionosphere  at  a height  above  ground,  A^  ► l/H^  and  A -*  1/2H  . 
Even  for  diffuse  ionospheric  boundaries  (as  treated  In  this  report),  the 
magnitudes  of  the  excitation  factors  are  of  the  same  order  as  the  reci- 
procal of  the  nominal  ionospheric  reflection  heights.  This  definition 
differs  from  that  often  found  in  the  literature,  where  the  nominal 
reflection  height  is  factored  out  and  appears  explicitly  in  equations 
analogous  to  F.qs.  (1)  and  (2),  and  the  excitation  factors  are  treated 
as  dimensionless  quantities  with  order-of-magnltude  unity  rather  than 
1/H  . Stated  differently,  in  F.q.  (1)  corresponds  very  roughly  to 
the  ratio  A /ll  often  used  in  treatments  that  model  the  earth- ionosphere 
waveguide  as  having  a well-defined  reflection  height. 


NUMERICAL  RESULTS 

All  results  given  below  were  calculated  using  PSR's  full-wave 
propagation  code  described  by  Field,  et  nl.  (19? 6*).  Table  1 gives  the 
attenuation  rate,  n,  relative  phase  velocity,  (v/c)  -1,  and  excitation 
factor.  A,  for  frequencies  of  45  Hz,  75  Hz,  20  kHz,  and  35  kHz,  and 
the  five  daytime  model  ionospheres  described  above.  A ground  conduc- 
tivity of  4 mhos/m  is  assumed  in  all  cases.  Results  for  only  the 
lowest  mode  (No.  1)  are  shown  for  45  Hz  and  75  Hz,  because  higher-order 
ELF  modes  are  well  below  the  cutoff  frequency  of  the  ear t h- ionosphere 
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Table  1 


CALCULATED  ELF/VLF  MODE  PARAMETERS 


Freq.  (Hz) 


Mode  No. 


Model 


r?  (dB/Mm) 


45 

1 

Standard 

Bain 

A 

B 

* 

C 

-7.99x10  1 
-7.29x10  1 
-7.27x10  1 
-6.61x10  1 
-6. 86x10'’ 

1.86x10  1 
' 2.18x10’’ 
-2.03  x 10' 1 
•2.13x10  1 
-1.73x10'’ 

9.23x10® 
9.95x  10  6 
9.62x10  6 
9.84  xlO'6 
8.93x10'6 

Standard 

-1.14 

-1.68x10'’ 

9.00  xlO'6 

Bain 

-1.15 

-2.03x10  1 

9 78x10 6 

75 

1 

A 

-1.07 

-1.88x10  1 

9.46x10  6 

B 

•1.04 

-1.99x10  1 

9.72x106 

C 

9.34  x 10' 1 

-1.60x10'’ 

8.82x10® 

Standard 

-1.88 

1.64x10  3 

1.07x10'® 

Bain 

3.26 

1.49x10  3 

1.13x10  5 

1 

A 

-2.44 

-1.41  x 10'3 

1.15x10® 

B 

2.78 

•1.38x10  3 

1.16x10  5 

C 

-2.00 

-1.55x10  3 

1.10x10® 

20  * 103 

Standard 

-7.97 

1.07x10  2 

1.78x10® 

Bain 

•14.71 

1.06X10'2 

1.71x10'® 

2 

A 

•11.29 

1.16X10'2 

1.79x10  5 

B 

•12.77 

1.15X102 

1.72x10  5 

C 

8.83 

1.1 1 x 10'2 

1.80x  10  5 

Standard 

-3.32 

-3.54  xlO'3 

2.99x10® 

Bain 

-6.39 

3.02x10  3 

4.99x10® 

1 

A 

4.35 

-3.14x10  3 

4.28x  10  6 

B 

5.31 

-3.08x  10'3 

4.60x10® 

C 

-3.47 

3.39x10  3 

3.42x10® 

J J * IU 

Standard 

4.53 

8.49x  10  4 

1.97x10® 

Bain 

•12.05 

1 . 1 3 x 10  3 

2.07x10® 

2 

A 

7.31 

1.32X10'3 

2.08x10  ® 

B 

9.10 

1.25x10  3 

2.05x10  5 

C 

-5.14 

1.05x10  3 

2.0  2x  10'® 

v/c  - 1 
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cavity.  For  VLF  (20  kHz)  and  LF  (35  kHz),  results  are  shown  for  the  first 
two  TM-waveguide  modes.  At  higher  frequencies,  the  attenuation  rates  of 
the  various  modes  are  nearly  equal,  with  the  result  that  a prohibitive 
number  of  modes  must  be  retained  to  achieve  reasonable  accuracy.  At 
these  higher  frequencies,  a geometric-optical  approach  is  generally  more 
practical  than  the  normal-mode  approach  used  here. 

The  results  shown  In  Table  l are  used  in  conjunction  with  Eqs.  (1) 
and  (2)  to  calculate  the  field-strength  versus  distance  graphs  shown 
in  Figs.  3 through  b.  Distances  shorter  than  several  hundred  mega- 
meters are  not  shown  on  the  graphs,  because  the  number  of  modes  used 
inadequately  represents  short-path  propagation.  Even  the  results  shown 
do  not  achieve  high  accuracy  until  the  path  length  exceeds,  say,  one  or 
two  megameters.  Because  we  are  interested  solely  in  the  relative  values 
of  the  field  for  the  various  models,  an  arbitrary  amplitude  scale  is 
used  in  Figs.  3 through  6. 

To  assist  in  the  interpretation  of  the  results  shown  in  Figs.  3 
through  6,  we  have  calculated  the  normal  reflection  coefficient,  mRm, 
versus  the  cosine,  C,  of  the  angle  of  incidence  for  each  of  the  five 
ionospheric  models.  Because  the  concept  of  a reflection  coefficient  is 
not  useful  at  F.LF,  where  the  wavelength  greatly  exceeds  the  distance 
between  the  earth  and  ionosphere,  |||R|||  is  given  only  for  frequencies 
of  20  kHz  and  35  kHz  (Figs.  7 and  8). 

As  an  additional  diagnostic  tool.  Figs.  9 and  10  show  the  relative 
Joule-heating  height-profile  for  the  lowest  waveguide  mode  at  75  Hz  and 
20  kHz.  The  figures  compare  the  heating  profiles  between  the  standard 
ionosphere  and  the  ionosphere  having  the  "Bain"  C-layer  ledge.  Note 
that  the  shape,  rather  than  the  magnitudes,  of  the  heating  profiles 
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should  be  compared  because  all  profiles  have  been  normalized  to  have  a 
peak  value  of  unity.  The  heating  profiles  above  about  70  km  are  less 
accurate  than  at  lower  altitudes  because  of  the  omission  of  the  geomagnetic 
field  in  the  calculations. 

DISCUSSION 

The  results  in  Table  1 for  45  Hz  and  75  Hz  show  that  the  effects  of 
the  C-layer  ledge  on  ELF  propagation  are  small.  Only  slight  differences 
in  attenuation  rate,  n,  and  excitation  factor.  A,  are  noted  among  the 
different  models  considered.  Moreover,  no  clear  trend  is  evident — 

Model  B giving  the  lowest  attenuation  at  45  Hz  and  Model  C the  lowest 
attenuation  at  75  Hz.  Figures  3 and  4 give  the  ELF  field-strength  versus 
distances  for  these  frequencies,  and  again  indicate  only  minute  differ- 
ences for  the  various  ionospheric  models  used.  Even  at  a propagation 
distance  of  10  the  spread  in  calculated  field  strength  is  only  about 
35  percent  (2  dB) . 

The  results  for  VLF/LF  are  very  different  than  those  described 
above  at  ELF.  Table  1 shows  that  at  20  kHz,  significant  differences  in 
calculated  attenuation  rates  occur  among  the  various  models.  Even 
larger  differences  in  n are  shown  at  LF  (35  kHz).  However,  no  important 
dependence  of  the  excitation  factor  on  the  C-layer  ledge  is  noted. 

The  trend  in  the  results  for  VLF/LF  is  much  more  regular  than  at 
ELF.  In  every  instance,  the  standard  DNA  model  gives  the  lowest  attenua- 
tion rate,  indicating  that  the  effect  of  the  C-layer  ledge  assumed  here 
is  to  degrade  the  propagation.  For  all  cases,  the  Bain  model  of  the 
C-layer  ledge  gives  the  highest  attenuation  rate. 


e 


1 


The  results  for  Models  A,  B,  and  C are  of  Interest  because  they 
Illustrate  the  dependence  of  propagation  on  the  altitude  of  the  layer — 
the  electron  density  and  thickness  having  been  held  constant  for  these 
three  models.  For  all  cases,  the  lowest-altltude  model  (B)  gives  the 
highest  attenuation  whereas  the  highest-altitude  model  (C)  gives  attenua- 
tion rates  only  slightly  greater  than  does  the  standard  model.  Two 
competing  effects  occur  when  the  altitude  of  the  layer  is  changed.  First, 
because  of  the  decrease  of  collision  frequency  with  altitude,  a given 
electron  density  corresponds  to  a greater  refractive  index  at  say, 
bl  km  than  at  55  km.  Second,  and  conversely,  a 100  el/cm^  increase  in 
electron  density  is  a relatively  modest  addition  to  the  standard  density 
at  61  km,  but  is  a very  large  relative  enhancement  at  55  km.  For  the 
models  assumed,  the  second  effect  is  clearly  dominant  because  the  ledge 
centered  at  55  km  caused  the  greatest  increase  in  attenuation  rate.  Of 
course,  this  trend  of  the  effect  of  the  ledge  increasing  as  its  altitude 
is  decreased  would  not  continue  indefinitely.  Placing  the  ledge  at  too 
low  an  altitude  (e.g.,  40  km)  would  have  virtually  no  effect  on  propaga- 
tion, because  the  electron  mobility  would  be  destroyed  by  the  large 
collision  frequency. 

Figures  5 and  6 show  that  the  effects  of  a C-layer  ledge  on  VLF/LF 
field  strength  can  be  considerable,  provided  that  the  ledge  exists  over 
lateral  distances  of  many  megameters.  At  20  kHz,  Fig.  5 shows  that  the 
spread  among  the  calculated  fields  is  about  a factor  of  four,  which 
corresponds  to  12  dB.  At  35  kHz,  Fig.  6 shows  spreads  of  up  to  a factor 
of  40,  which  corresponds  to  about  30  dB. 

Figures  7 and  8 provide  insight  into  the  behavior  of  the  waveguide 
modes  given  above.  Note  that  the  full-wave  mode  calculations  indicate 
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that  the  eigen-cosine  of  the  lowest  VLF/hF  modes  are  about  0.2  at  20  kHz 
and  0.15  at  35  kHz.  These  eigen-cosines  correspond  physically  to  the 
cosine  of  the  inrid  ice  angle  at  the  ionosphere.  For  C " 0.2,  Fig.  7 
shows  that  the  ordering  of  | i |R| | | for  the  various  models  at  20  kHz  is 
the  same  as  the  ordering  of  n in  Table  l;  i.e.,  the  reflection  coeffi- 
cient Is  lowest  (n  highest)  for  the  Bain  model,  and  so  forth.  Similar 
conclusions  are  drawn  from  Fig.  8 for  35  kHz  and  C * 0.15. 

An  Important  conclusion  drawn  from  Figs.  7 and  8 is  that  the  relative 
effects  of  the  various  assumed  C-layer  ledges  depend  strongly  on  inci- 
dence angle.  For  example,  as  discussed  above,  the  Bain  model  gives  the 
lowest  reflection  coefficient  for  long  propagation  paths,  which  corres- 
pond at  oblique  incidence  (C  £ 0.2).  For  steeper  incidence  angles  (C  ^ 0.6; 

0 - 53°),  however,  the  Bain  model  gives  the  highest  reflection  coefficient  of 
any  of  the  models  considered. 

Figure  9 indicates  why  the  C-layer  ledge  has  such  a small  influence 
on  ELF  propagation.  As  shown  by  the  Joule-heating  profile,  the  important 
ionospheric  Interactions  occur  over  a broad  range  of  altitudes  that 
extend  from  about  50  km  to  85  km  at  75  kHz.  Only  the  lowest  10  km  or 
so  of  this  range  is  significantly  affected  by  the  Bain  model  of  the 
C-layer  ledge.  Moreover,  the  enhanced  heating  between  50  km  to  55  km 
Is  largely  compensated  by  reduced  heating  between  55  km  to  60  km, 
causing  the  attenuation  rate — which  is  essentially  proportional  to  the 
area  under  the  heating  profile — to  he  relatively  unaffected. 

Conversely,  at  20  kHz,  Fig.  10  shows  that  for  the  standard  model  the 
heating  peak  occurs  in  a well-defined  region  about  1 5-km  thick,  centered 
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at  about  bS  km.  The  main  effect  of  the  Bain  C-layer  Is  to  lower  this 
heating  region — and,  presumably,  the  reflection  height — by  nearly  10  km. 
Although  not  evident  from  the  normalised  curves  shown  In  Ktg.  10,  the 
heating  peak — and  the  area  under  the  heating  curve — for  the  Bain  model 
Is  about  twice  as  large  as  for  the  standard  model. 
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III.  RELATIONSHIPS  AMONG  INCIDENCE  ANGLE , 

FREQUENCY,  AND  VLF/LF  REFT  .KCTTON  HE TCI  ITS 

The  range  of  altitudes  from  which  VLF/LF  signals  are  reflected  from 
the  ionosphere  depends  on  the  incidence  angle  and  frequency  of  the  wave. 
Ordinarily,  the  reflection  height-range  tends  to  be  lowered  if  1)  the 
wave  frequency  is  reduced,  or  2)  if  the  incidence  angle  is  made  more 
oblique.  Thus,  for  a given  frequency,  long-path  signals  tend  to  interact 
with  lower  ionospheric  heights  than  do  short-path  signals,  which  have 
steep  incidence  angles.  For  the  reasons  stated  above,  however,  it  should 
be  possible  to  cause  oblique  waves  to  interact  with  the  same  ionospheric 
regions  as  steep-incidence  waves  by  properly  adjusting  the  frequency. 
Accordingly,  this  section  analyzes  the  relationship  between  incidence 
angle  and  frequency  of  VLF/LF  waves  that  reflect  from  a given  ionospheric 
height-range. 

BACKGROUND 

Some  progress  on  the  above  problem  has  already  been  reported  in  the 
literature.  Booker  and  Crain  (1964)  computed  reflection  loss  by  assuming 
(on  intuitive  grounds  without  proof)  that  VLF/LF  reflection  takes  place 
near  a level  where 
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2uicos  0 **  to  /v  5 u)„ 
p e 0 


(3) 


In  F.q.  (3),  w is  the  angular  frequency  of  the  wave,  0 is  the  incidence 

2 

angle  at  the  ionosphere,  vg  is  the  electron-collision  frequency,  and 
is  the  electron-plasma  frequency,  which  depends  linearly  on  the  electron 
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density,  N^.  Thus,  the  characteristic  frequency,  depends  on  height 

through  the  ratio,  N /v  ; and  Eq.  (3)  states  that  two  waves  will  sample 

6 6 

the  same  height-range  (i.e.,  the  same  values  of  N^/v^)  provided  that  the 

2 

product  iocos  0 is  kept  constant. 

Because  Booker  and  Crain  simply  wrote  down — rather  than  derived — Eq. 
(3),  they  give  no  criteria  for  its  validity  or  range  of  applicability. 

Of  course,  given  the  absence  of  the  gyrofrequency  in  Eq.  (3),  the  applica- 
bility is  clearly  limited  to  daytime  or  disturbed  conditions  where  geo- 
magnetic effects  play  a relatively  minor  role.  Equation  (3)  was  put  on 
a more  rigorous  footing  by  Field  and  Engle  (1965),  who  show  that  provided 


and 


cos  0 <<  1 


w 0OCe 


BZ 


(4a) 


(4b) 


the  main  reflections  occur  in  a several-kilometer  altitude  range  centered 
near  the  altitude  where 


2 

/2  tocos  0 » io^  , (5) 

which  is  essentially  the  same  as  Eq.  (3),  and — like  Eq.  (3) — is  invarient 

2 

to  changes  in  to  or  cobO  provided  the  product  tocos  0 remains  constant. 
DERIVATIONS 

The  task  at  hand  is  to  derive  relations  analogous  to  Eqs.  (3)  or  (5) 
that  are  not  restricted  by  the  conditions  (4a,  4b)  and  that  apply  under 
conditions  in  which  the  geomagnetic  field  must  be  included.  The  approach 
taken  is  to  derive  expressions  for  the  ionospheric  reflectivitv-per-unit- 
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height,  R(z) , and  determine  which,  if  any,  combinations  of  id  and  cosO 
will  leave  R unchanged.  Only  partial  success  is  achieved  because,  a] though 
a useful  relation  unrestricted  by  conditions  4a,  4b  is  found,  the  equa- 
tions for  conditions  where  the  geomagnetic  field  must  be  included  are  too 
complicated  to  oe  of  practical  value. 


Vertical  Polariz atlon 

We  begin  by  defining  the  following  three  quantities  that  occur  through- 
out the  derivation: 


U) 

2 . A 4 0 
n as  l + i — 

(A) 


2 2a  » < 0 

q s cos  0 + i — 


Q - q/n 


(b) 

(7) 


In  the  above  expressions,  ui^,  which  is  a function  of  height,  is  defined  by 
2 

Kq.  (3),  n is  the  ionospheric  refractive  index,  and  q is  the  solution  to 
the  well-known  Booker  quartic.  Equations  (6)  and  (7)  are  based  on  the 
approximation  that 


v > u>  and  u>  , 

e c 

which  is  valid  below  70  km  or  so  for  VLF/LF. 

The  exact,  full-wave  equation  for  the  reflection  coefficient  is  non- 
linear and,  hence,  must  be  solved  numerically.  Because  an  analytic 
relation  between  i»>  and  cos9  is  desired,  some  approximations  must  be  made. 
Specifically,  attention  is  restricted  to  conditions  where  the  reflectivity- 
per-un it -height  is  not  too  large. 

For  vertical  polarization,  the  relevant  Maxwell  equations  are 
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~ = -ikqQH  (9) 

and 

H = -ikn2E  , (10) 

where  k is  the  free-space  number,  E is  the  horizontal  electric  field,  and 
H is  the  normalized  horizontal  magnetic  intensity  ( Budden , 1961 , p.  142). 
The  method  of  solution  is  based  upon  deriving  correction  terms  to  the  WKB 
approximation,  which  neglects  reflections  altogether.  Thus,  we  insert 


u,d 


(11) 


into  Eqs.  (9)  and  (10),  where  the  subscripts  u,d  denote  upgoing  and  down- 
going waves,  respectively.  After  some  rearrangement,  the  following  coupled 
equations  for  the  upgoing  and  downgoing  waves  are  obtained: 

h;+[^+Hhu'Shi  <12a) 

Hd  +[2Q~lk<1]  Hd  " 2QHu  ’ ^12b^ 

where  the  prime  denotes  z-dif ferentiation.  The  quantity  Q'/2Q  represents 
coupling  between  upgoing  and  downgoing  waves,  and  is  treated  as  a perturba- 
tion term. 

The  zero-order  solution  to  Eq.  (12a),  obtained  by  neglecting  the  right- 
hand  side,  is  easily  shown  to  be 


Specifically,  H differs  from  the  unnormalized  magnetic  intensity  by 
a factor  of  12011  ohms. 


H CO80 

H - exp 

U /Q 
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q(C)dC 


where  is  the  wave  Incident  on  the  ionosphere  from  below;  i.e.t 

at  z « 0.  Equation  (13)  is  the  well-known  WKB  representation  of  an  upgoing, 

vertically  polarized  wave. 

The  downgoing  wave  generated  by  the  upgoing  wave  is  calculated  to  first 
order  by  inserting  Eq.  (13)  into  the  right-hand  side  of  Eq.  (12b): 


Hd  + [2^  " lkqj  Hd  " 2^  exp|-ik  /q(°di:j 


. (14) 


Equation  (14)  is  a linear,  first-order,  differential  equation  that  may  be 
integrated  immediately.  The  resulting  expression  for  H^,  subject  to  the 
physically  necessary  condition  that  ■ 0 at  z ■ 00 , is 
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| 2q(0 


The  reflection  coefficient  at  the  ground,  Rv,  is  obtained  by  setting 
z ■ 0 in  Eq.  (15),  and  noting  that  Hd(0)  - RH^,  and  Q(0)  - cosO,  whence 


f r i c 

- J dC  exp  -2ikJ* qdc  , 

0 L L 0 
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which  is  a well-known  result  for  weakly  reflecting  media  derived  in  a 
different  manner  by  Wait  (1970). 

The  terms  in  Eq.  (15)  are  grouped  to  permit  a clear  physical  interpre- 
tation. The  first  curly  brackets  contain  the  expression  for  the  upgoing 
wave  at  a height,  c;  the  second  curly  brackets,  the  expression  for  the 
reduction  in  strength  of  a downgoing  wave  in  propagating  from  a height  f, 
to  a lower  height,  2.  Thus,  the  quantity  in  the  third  bracket  denotes  the 
amount  of  upgoing  wave  converted  into  downgoing  wave  in  the  thin  layer 
extending  from  t to  r,  + dc.  A similar  interpretation  can  be  obtained  from 
Eq.  (16)  by  noting  that  the  exponential  factor  in  the  integrand  accounts 
for  the  loss  suffered  by  an  upgoing  wave  propagating  to  the  height,  C;  and 
the  loss  suffered  by  a downgoing  wave  propagating  from  the  height,  C,  back 
down  to  the  ground  at  t - 0.  Thus,  the  reflectivity-per-unit-height , R^, 
can  be  written,  by  using  Eqs.  (6)  through  (8), 


2(tocos  6 


which  has  the  following  magnitude: 
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Horizontal  Polarization 

The  derivation  for  horizontal  polarization  proceeds  exactly  as  de- 
scribed above  for  vertical  polarization,  except  that  the  electric  field 
is  used  rather  than  H.  Without  going  into  details,  note  that,  aside  from 
unimportant  multiplicative  factors,  the  equations  for  horizontal  polari- 
zation are  Identical  to  those  for  vertical  polarization,  except  that  the 
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quantity,  Q,  Is  replaced  by  the  simpler  quantity,  q.  The  reflection 
coefficient  is 


(19) 


and  the  ref lect ivity-per-unit  height  is  given  by 


[fi] 


2^wcos“fl+  iu>pj 


which  has  the  magnitude 
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u>  cos  0 + u>q 

r 

(20) 


(21) 


DISCUSSION 

Subject  to  the  approximations  made  in  the  derivations,  and  R^ 
determine  the  heights  at  which  reflection  occurs.  Of  course,  as  indicated 
by  Eqs.  (16)  and  (19),  Ry  and  R^  do  not  completely  determine  the  ground- 
level  reflection  coefficients,  which  also  depend  on  the  exponential  term 
accounting  for  losses  suffered  propagating  to  and  from  the  reflection 
regions.  Thus,  the  following  subtle  distinction  arises.  If  one  wishes 
to  obtain  the  same  reflection  coefficient  for  obliquely  and  steeply  inci- 
dent waves,  then  w and  cos6  must  be  adjusted  such  that  the  entire  integrand 
of  Eqs.  (16)  or  (19)  remains  constant.  However,  if  the  goal  is  to  use  an 
obliquely  and  a steeply  incident  wave  to  sample  the  same  ionospheric  height 
ranges,  then  id  and  cosO  must  be  adjusted  in  a way  that  Ry  and  R^  remain 
unchanged.  In  the  latter  case,  the  ground-level  coefficients  will  not 


necessarily  be  the  same. 
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In  principle,  Rv  and  R^  account  for  both  gradient  reflection,  charac- 
terized by  large  values  of  u>^,  and  classical-turning-point  reflection, 
character iz led  by  either  q or  Q ^5 0.  Because  of  the  complex  functional 
forms  used  to  approximate  q or  Q (see  Eqs.  (7)  and  (8)),  neither  q or  Q 
can  vanish  in  the  altitude  range  considered  here,  and  only  gradient  reflec- 
tion is  a factor. 

In  Eqs.  (17),  (18),  (20),  and  (21),  the  height  dependences  of  and 

R^  are  contained  solely  in  ui^,  and  its  derivative,  u>^,  whereas  the  dependence 

on  (0  and  cosO  is  shown  explicitly.  The  question  to  be  answered  is  whether  a 

simple  functional  relationship  can  be  found  whereby  u>  and  cosO  can  be  varied 

without  altering  Rv  or  R^.  The  answer  is  remarkably  simple  for  horizontal 

polarization  because,  as  shown  by  Eq.  (20),  depends  on  frequency  and 

incidence  angle  solely  through  the  combination  ux?os“0.  Thus,  for  horizontal 

polarization,  the  conclusion  that  the  ref lect ivity-per-unit-height  is 

2 

invar ient  to  changes  in  <u  or  cosfl,  provided  that  weos  P remains  constant, 
applies  even  when  the  conditions  (4a)  and  (4b)  are  violated.  Of  course, 
this  conclusion  is  not  totally  unqualified,  because  the  derivation  is  based 
upon  the  assumption  that  the  medium  changes  slowly  enough  that  the  product 
of  R^  and  the  local  inverse  vertical  wave  number,  (kq)  is  less  than  unity; 
l.e.,  the  condition 


(22) 


must  hold.  Of  course,  Eq.  (22)  is  the  validity  criterion  for  the  WKB 

i 

approximation,  and  Field  and  Engle  (1964)  derived  Eq.  (5)  simply  by  maxi- 


mizing the  left-hand  side  of  Eq.  (22). 
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The  situation  is  more  complicated  for  vertical  polarization  because, 

as  shown  by  F.q.  (17),  depends  on  ui  and  cosO  through  factors  other  than 
> 

tocos*" 0 • For  either  oblique  or  steep  incidence,  however,  approximate  forms 
for  result,  from  which  relatively  simple  conclusions  can  be  drawn.  For 
sufficiently  oblique  angles,  F.q.  (17)  becomes 


*> 

cos“0  <<  S 


(23) 


whereas,  for  sufficiently  steep  angles,  Eq.  (17)  becomes 

s 

R % y if  sin“P  <<  S , (24) 

2(iocosZP  + 1u>q) 

which  is  identical  with  Eq.  (20)  because  the  distinction  between  vertical 
and  horizontal  polarization  disappears  for  nearly  normal  incidence.  Both 
Eq.  (23)  and  Eq.  (24)  depend  only  on  the  combination  tocos  6.  Further  note 

that  the  magnitude  of  R^  for  oblique  and  steep  incidence,  as  described  by 

* 

F.qs.  (23)  and  (24).  is  Identical  although  the  phases  are  different. 

We  thus  conclude  that  steeply  and  obliquely  incident,  vertically  polar- 
ized signals  are  reflected  from  about  the  same  ionospheric  heights  provided 

that  uicos‘0  is  held  constant.  This  conclusion  should  be  reasonably  accurate, 
2 2 

provided  cos  0 <<  S or  sin  0<<  S;  i.e.,  for  incidence  angles  greater  than 
about  70°  or  less  than  about  20°.  For  Intermediate  angles  and  vertical 
polarization,  no  such  relationship  exists,  as  can  be  seen  from  F.q.  (18). 


This  conclusion  also  follows  from  Eq. 


(18),  with  0*0  or  0*u/2. 
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This  section  has  been  concerned  with  ionospheric  conductivity  height- 
profiles  essentially  devoid  of  rapid  transition  regions,  which,  because  of 
large  height-derivatives,  cause  validity  criteria  such  as  F.q.  (22)  to  be 
violated.  It  is  precisely  at  such  "ledges,"  however,  that  waves  are  reflect- 
ed regardless  of  frequency  and  incidence  angle  (aside  from  the  Brewster 
angle).  Thus,  in  the  approximation  opposite  to  that  used  here — the  nearly 
sharply  hounded  layer — both  steeply  and  obliquely  incident  waves  will  reflect 
from  the  same  height  regardless  of  frequency,  provided  only  that  the  wave 
can  penetrate  to  and  from  the  layer. 

We  were  not  able  to  derive  useful  relations  between  to  and  cost1  for 
situations  in  which  the  geomagnetic  field  must  be  included.  Aside  from 
very  special  conditions  in  which  the  propagation  is  essentially  isotropic 
(e.g.,  east-west  propagation  at  the  equator,  normally  incident  circularly 
polarized  waves  at  the  geomagnetic  poles),  the  equation  for  reflectivity- 
per-unit-height  becomes  much  too  complicated  to  be  of  practical  analytic 
value . 
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